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Abstract 

We establish a novel representation of arbitrary Euler-Zagier 
sums in terms of weighted vacuum graphs. This representation 
uses a toy quantum field theory with infinitely many propagators 
and interaction vertices. The propagators involve Bernoulli poly- 
nomials and Clausen functions to arbitrary orders. The Feynman 
integrals of this model can be decomposed in terms of a vertex 
algebra whose structure we investigate. We derive a large class of 
relations between multiple zeta values, of arbitrary lengths and 
weights, using only a certain set of graphical manipulations on 
Feynman diagrams. Further uses and possible generalizations of 
the model are pointed out. 



1 Introduction 



The perturbative evaluation of Green's functions in quantum field theory 
lead to a class of iterated parameter integrals whose explicit calculation 
becomes very difficult beyond the first few orders in the coupling constant 
expansion. Any progress in this area of work must be based on an intimate 
knowledge of the properties of various types of special functions such as 
polylogarithms, hypergeometric functions, and their generalizations (see e.g. 



In recent years, some structure is seen to emerge from the seemingly 
haphazard occurence of those special functions as the values of individual 
Feynman diagrams. Kreimer's hypothesis based on a rule of associating 
knots to Feynman diagrams, allows one to predict from knot-theoretical 
considerations the level of transcendentality which can possibly appear in the 
counterterm coefficients of an ultraviolet divergent diagram. Even though 



it has been verified for a large number of examples [13| the raison d'etre for 
the correspondence between graphs and knots remains presently mysterious. 
More recently there are indications that knot-theoretical concepts may be 



of relevance even for the finite parts of Feynman diagrams [11|. 

The remarkably rich mathematical structures surfacing in this corre- 
spondence make Feynman diagrams increasingly interesting from the pure 
mathematician's point of view. The objects encountered in the calculation 
of UV divergences in perturbative quantum field theory, multiple harmonic 
sums, are of considerable relevance to number theory and other branches of 



mathematics (see e.g. |5|, ||, g |l9|, |g, |T|, ||] ) 



Quantum field theory amplitudes can be calculated in coordinate space 
or in momentum space. In four-dimensional field theory the arising integrals 
are normally of a similar type and degree of difficulty. This is very different in 
the case of a one-dimensional quantum field theory compactified on a circle, 
which will be considered in the present paper. Such quantum field theories 
arise naturally if one represents one-loop amplitudes in iD-dimensional field 
theory in terms of first-quantized path integrals. An approach to quantum 
field theory along these lines has gained some popularity in recent years 
after it was discovered that it allows one to reorganise ordinary field theory 



amplitudes in a manner similar to string theory amplitudes p2\ , ^ |3^, 36, ^ . 
In this type of formalism the D-dimensional space-time enters as a target 
space, and amplitudes are calculated in terms of an auxiliary field theory in 
one-dimensional parameter space. Green's functions in parameter space are 
then used for the evaluation of Feynman diagrams in this one-dimensional 
'worldloop' theory. 
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As a simple example, let us consider the one-loop effective action for 
a scalar field theory with a interaction. This effective action can be 
expressed in terms of a first-quantized path integral as follows, 



rw = ir^e--/ ^.Me-^" (1.1) 

^Jo J- Jx{T)=x{0) 

Here T is the usual Schwinger proper-time for the particle circulating in 
the loop. At fixed T a path integral has to be calculated over the space of 
closed loops in spacetime with period T. This integral contains a zero mode 
which is removed by fixing the center-of-mass of the loop = y drx(r). 
The reduced path integral is evaluated perturbatively by expanding the 
interaction exponential and using the parameter space Green's function |^ 



Momentum space methods have also sometimes been used, which in this 
case lead to Fourier sum representations. In |3|, [l^ a number of such sums 
were calculated to provide a check on the resolution of certain ambiguities 
which in the first-quantized formalism can arise in curved backgrounds. In 
this comparison one finds that terms given by simple polynomial integrals 
in coordinate space may, in momentum space, correspond to non-trivial 
multiple sums of the Euler-Zagier type. 

In the present work we will turn the logic around, and use this formalism 
as a tool for the systematic study of Euler-Zagier sums. Euler-Zagier sums, 
also called multiple C, values or multiple harmonic series, are defined by 

C{ki,...,kra) = ^ fc^ ^ =Lifc^...fc^(l,...,l). 



(1.3) 



They are special values of the multidimensional polylogarithms Li^.^^ fc^, 
defined as 

Lifc,..,^(zi,...,z^) ^ X X • (1-4) 

ni>n2> - >n,^>0 "'I "' ""^ 



'^The constant part of this Green's function is irrelevant for the final physical results 
and usually deleted from the beginning. 
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One calls m the length (or depth) of such a series, and ki + k2 + • • • + km 
its level (or weight). Sums of the type ( |1.3D were first considered by Euler 
[p!6| . Euler himself noted that numerous relations exist between Euler-Zagier 
sums. Some simple examples are the following (all given by Euler) 

C(2,l) = C(3), (1.5) 
C(3,l) = ^C(4)-^C^(2) = ^, (1.6) 

C(2,2) = ^C^(2) 4C(4) = ^, (1.7) 

C(3,2) = -^c(5) + 3C(2)C(3), (1.8) 
C(4,l) = 2C(5)-C(2)C(3). (1.9) 

Further results for the length two case can be found in ||4^, |l|. Systematic 
investigations of Euler-Zagier sums of length higher than two have been 
undertaken only in recent years ||, ||, |o|, |2|, H, |, |2|, |, [S^, |, H. 



From the point of view of physics the study of their relations is relevant for 
attempts at a classification of the possible ultraviolet divergences in quantum 
field theory plllll. 



We would like to be able to represent arbitrary such sums in terms of one- 
dimensional Feynman diagrams. To achieve this goal we have to generalize 
the usual worldline path integral formalism in the following ways: 

1. As explained above, the first-quantized loop path integral is defined to 
run over the space of all periodic functions, with the constant functions 
eliminated. Here we will restrict it to one "chiral half" spanned by 
the basis functions fn{u) = e^'^'™,n = 1,2,.... (This amounts to a 
complexification of spacetime.) 

2. We choose the kinetic term of our model in such a way that arbitrary 
inverse powers of derivatives will appear. 

Those purely mathematical considerations lead us to define the "C-model", 
a one-dimensional quantum field theory given by the following partition 
function, 

Zig,X) = [ Vx{u)e~^, 
Jh 

S= [ dui [ du2lx{ui)(l-X2md-Ax{u2)- [ du e5^(")+^^(") . 

(1.10) 
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Here and in the following d = denotes the ordinary derivative. The path 
integral is to be performed over the Hilbert space 

oo oo 

?^ = {x(n)|x(u) = ^a„e2™", ^ |a„p < oo}. (1.11) 

n=l 71=1 

The perturbative expansion of both the kinetic and the interaction terms 
for this toy model leads to the following Feynman rules, 



Vertices V^''^: 



Propagators: 1 




p,q = 0,l,2,... 



2 /c = 0,l,2,... 



Figure 1: Feynman rules of the C-model. 

(k) 

In section 2 we write the propagators explicitly in terms of Bernoulli 
polynomials and Clausen functions. We show that any multiple C sum ( |l.3| ) 
can be represented as a Feynman diagram in this model. In section 3 we 
investigate the properties of the elementary tree-level n-point integrals with 
arbitrary external propagators. Section 4 demonstrates how one can use 
partial integrations and reality conditions to derive a large class of relations 
between multiple ( sums. In section 5 we point out possible further uses of 
the model, as well as some generalizations. 



2 Basic Properties 

Inverting the kinetic part of our (non-local) Lagrangian ( l.lOp , and writing 
in the defining basis of the Hilbert space Ti, we find 

. g^'Hu,.) = E--k- (2-1) 

n=l (2™) 

If we represent the unit circle in the complex plane, this sum will turn into 
the k-th polylogarithm. 
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(ui2 = ui — U2, Zi = e^™^). We note the following properties of ^f^'^^ 
d (k) d (k) (fc-i) 

= {-lW=g^'\l-uu), (2.4) 
.^^)(0) ^ (2.5) 

/ dni,2 5i? = 0. (2.6) 
Jo 

(p. 3D can be inverted using the explicitly known integral kernel for inverse 



derivatives in this space |36], 



git"'^ = C du{ui\d-'\u)g^^\u-U2)^ (2.7) 
Jo 

(ni|a-"|u2) = -i-i?„(|ni2|)sign"(^.i2) (2.8) 
n! 

n-l 



Bn{ui2) , , . , N -.N „\ 

= + 77^^^ sign^i2 - 1 . 2.9 

n! l[n — \)\ 

(k) 

Here Bk denotes the fc-th Bernoulli polynomial. To write gj^2 more explicitly 
we split it into its real and imaginary parts, 

ffS? = (2.10) 
Using the standard integral representation of the polylogarithm 

it is then easy to show that 
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3n = ^('5(«i2)-l + icot(inii2)), (2.12) 



1 /I 

2 



(^^sign(ui2) - ni2 + ;^ln|2sin(7r-ui2)|^, (2.13) 



(2) l/i I 2 1 , i /"^ ln^sin(27rni2) 
5i2 = 7 1^121 -Ui2 - - + ^ / 



4V' 6 7r2yo ^ 1 -2ecos(27rni2) 



2.14) 

- 2U12)(|'U12| - W12) 

^ dC 



5'l2'' = -^(sign(ui2) - 2ui2)(|ui2| - n?2 



+ ylnein(l-2Ccos(27rni2)+^'), 

(2.15) 



(k^O even) 11,, , 

9i2 = ~2A?!^^'''^^' 



fc 



+ ^fe^rd^ln-k sin(2..,2) 



(2vr)M^-l)! Jo M-2ecos(27r«i2) + e2' 

(2.16) 



(fc^l odd) 1 1 D /I l\ ■ / \ 

9l2 = -^j^Bk{\ui2\)sign{ui2) 



(2.17) 



{ui2 G [—1,1]). Note also that the imaginary part of §['2 is, up to a nor- 
mahzation factor, identical with the k-th Clausen function Clfe(27rui2) (see, 
e.g., 1^^). Only the real parts are present in the calculation of worldline 
path integrals in standard field theory. In particular, note that g^"^^ is, up to 
a conventional factor of 4, identical with the function G introduced in ( |1.2| ) 
(for T = 1). Only the imaginary parts are capable of producing C('T')'s with 
n odd. 

A vacuum diagram in our model involves a multiple integral on the unit 
circle with an integrand which is a product of propagators ( |2.1| ). (One of 
the integrations is redundant due to translation invariance.) The result of 
the integrations can obviously be decomposed as a sum of multiple C values. 
Our model thus defines a map from the set of weighted vacuum graphs to the 
(vector space generated by) the Euler-Zagier sums. This map is, moreover, 
easily seen to be surjective; the reader will immediately convince herself 
that with the above Feynman rules the following "sea shell" diagram (fig. 
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^) evaluates to 

\>^'^iggf'^~'C{ki,...,km). (2.18) 

The u-integrations produce 5 functions for "momentum conservation" at 
every vertex, while the Fourier sums of the inserted g^^^ propagators yield 
Heaviside step functions leading to the desired ordering of the remaining 
sums. 




Figure 2: 'Sea shell' diagram representing the general Euler-Zagier sum. 



3 The elementary vertex integrals 

Let us now start on an investigation of the properties of the Feynman in- 
tegrals in x{= M)-space. An obvious first step is to consider the folding of 
the elementary vertices with arbitrary sets of propagators. We denote the 
elementary vertex integral by 




Up — u) 



xg^'^\u - Up+i) • • • g^'''\u - Up+g). (3.1) 
We note that it has the following obvious properties, 

= ifp = 0or<7 = 0, (3.2) 

itX = (3-3) 
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3.1 Two- vertex integral 

By construction two-point vertices can be integrated out trivially, 

/'d-3<7f^g=^^r^ (3.4) 
Jo 

3.2 Three- vertex integral 

The evaluation of vertex integrals is complicated by singularities which can 
appear due to the presence of the cotangent function in g^^^ . Integrals involv- 
ing cot(7rui2) need to be performed using the principal value prescription. 
One way of calculating them is to transform them into complex contour 
integrals via the substitution z = exp(27rin), 

/d„nco.w„-«))=i''/^ni±|. (3.5) 

fc=l k=l " 

Those are evaluated by means of residues (the poles on the contour give half 
values due to the principal value prescription) ^ The first few integrals are 

1 

J du cot{TT{u — ui)) = 0, (3-6) 



1 2 

J dti cot(7r(n — Ujfc)) = 5{ui — U2) — 1, (3-7) 







k=l 



/ JJ cot(7r(n — u,t)) = (5(ni — n2) cot(7r(n2 — Us)) + 

-I- (5(ui - U3) cot(7r(ni - ti2)) + 



k=l 



+ 6 {u2 - U3) cot {it {U2 - Ul)). (3.8) 

Alternatively one can also calculate those integrals recursively using, under 
the integral, the following identity, 

cot(7rui2) cot(7rni3) -|- cot(7rti2i) cot(7rn23) + 

-|-cot(7rti3i)cot(7rn32) = -1 + 6{ui2)S{uis). (3.9) 



^ We re mark that, if one expresses the propagators of the sea shell diagram via (|2.2[) 



and (2.11), and then calculates the Ui-integrals by means of residues, then one arrives 



precisely at Kontsevich's integral representation [M for the multiple sum. 



8 



This identity will be of further use later on. With these results the vertex 
integral of three propagators g^^^ is 

1 

Io°o(ni,n2,n3) = J du g(°\ui - u)g^''Hu2 - u)g(^Hu - u^) 


= ^[^{U13)S{U23) - S{ui3) - S{U23) + 1 

+ {S{ui3) - l)icot(7rn23) 
+{S{u23) - l)icot(7rni3) 

- cot(7rti23) cot(7rtii3)] , (3.10) 

and can be identified as 

Ioo{ui,U2,U3) = g'f^g^^s- (^-H) 
Applying the identity (|3.9| ) gives also 



iQ0\Ui,U2,U3) - 512 523 -512 513 ~ 2^13 

_ (0) (0) _ (0) (0) _ 1 (0) 

— 521 5l3 521 523 2^23 

- 5i2 523 + 52i5i3 -2^2^13- \'^-^^) 

Folding of eq. ( p. Ill ) with {uii\d^'^\ui) and (^2' |9~™|n2) leads to 

1 

/0^(ni,n2,n3) = [ du g^^'^m - u)g^"^'\u2 - u)g^''\u - U3) 



9i^95\ (3.13) 



whereas the forms ( 3.12 ) can be folded with {u3\d {uy) 

1 

I^o{uuU2,U3) = j dug^^\ui-u)g^^\u2-u)g^^\u-U3) 


MAk) (0)(k) 1 (fc) 



5l2 523 5i2 5i3 2^13 

(0) (fc) _ (0) (k) _ 1 (fc) 
521 5l3 521 523 2^23 

5l2 523 + 521 5l3 - 2 ^2^13 • (.3.14) 
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Nontrivial is the case 

4fe(ui,U2,U3) = / dug^^\u-ui)g^''^u-U2)g^^Hu3-u) 
Jo 

= Uik{z3i,Zi2). (3.15) 

The general case of ah three indices different from zero can be reduced to 
this case by partial integrations, e.g. 

\ugW(n4i)g^^\n42)g^'\us4) 











+ l\ug('\n4i)g^'Hn42)9^^Hus4) 







= U2l{z3i,Zi2) +U2l{z32,Z2l)- (3.16) 

3.3 A three — point relation 

Beyond the three-point case the systematic investigation of the elementary 
vertex integrals becomes cumbersome. We will be satisfied here to note 
that, from the representation ( |2.12| ) for g^^^ and the identity (3^) one easily 



derives the following pair of (complex conjugate) three-point identities, 

n(0)„(0) , „(0) (0) (0) (0) _ , . JO) r „{0) , r „(0) r r 
521 ffsi +512 532 +513 523 " 1 + <J12fi'32 +031521 + 0235l3 "012013, 

5l2 5l3 + 521 523 + 531 532 " + 012^23 + 03l5l2 + 023531 " Ol20l3- 



(3.17) 




We may represent the first identity graphically as in fig. ^. 
111 1=3 1 

+X + X =1+ 

/ ^ ^ \ ^ / \ • 

3 2 3 2 3 2=1 3 2 3=2 1=2=3 

Figure 3: Three-point relation. 



The second identity has all arrows reversed. Those identities can be used 
to transform any vertex integral involving two g^^^^s which are either both 
ingoing or both outgoing. 
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By iteration of this three-point identity one can contruct an analogous 
identity for an arbitrary number of points. The formulas are rather cum- 
bersome and will not be given here. 

4 Derivation of Multiple ( Relations 

We will now show how can use the formalism developed above for deriving 
a large class of multiple C, relations by simple manipulations on graphs, with 
no need to ever explicitly write down sums. In those manipulations we will 
make use of the following elements: 

1. The triviality of the real part of g^°\ gf^ = \{5i2 - 1) (see (^l2D). 

2. Partial integrations. 

3. The three-point relations eq. ( p.l7 ). 

4. The vanishing of diagrams containing a vertex with only ingoing or 
only outgoing propagators. 

5. The two- vertex integration formula eq. ( |3.4| ). 

At intermediate steps sometimes ill-defined sums will appear such as 
Those will always cancel out in the final results. The appearance of 
divergent sums could be easily avoided using a regularization such as in |^] 
but we will not do so here. 

4.1 Length two 

Let us begin at length two. The simplest way of arriving at identities is 
provided by the first one of the points listed above. Consider the sea shell 
diagram at length two (fig. |^), representing C(a,6), and the same diagram 
with the middle propagator reversed (fig. representing ("(6, a). (We will 
generally disregard the coupling constant factors in the following.) 
Adding up both diagrams we can replace the middle propagator by twice its 
real part (fig. ^). Since = ^(^12 — 1) this diagram can then be replaced 
by the sum of the two diagrams shown in the right hand side of fig. ^. Using 
eq. ( ^.4| ) on the rightmost one of those we obtain the identity 

C(a,6)+C(6,a) = C(a)C(fe) - C(a + fe)- (4.1) 
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(a) 



(b) 



Figure 4: Diagrammatic representation of the reflection identity. 



This identity is weh-known ||2l|, |9[, and has been named "reflection formula" 
ini. 

Another way of obtaining identities is partial integration. Instead of 
considering the sea shell diagram (fig. ^a), which represents CiO'^b), let us 
consider the more general diagram (fig. |5|) that represents a number Gafi,c = 
Ga,c,b- For a = 0, 6 = or c = 0, it can be represented by multiple C 




Figure 5: Auxiliary diagrams for multiple ^ relations of length 2. 
functions, 

Go,b,c = mCic), Ga,b,o = Ga,o,b = C{a, b). (4.2) 

To see the first of these identities one adds the same diagram with the left 
propagator reversed, and uses elements one and four from the above list. For 
the later discussion we refer to these diagrams as zeta diagrams in contrast 
to the other, non-zeta diagrams. If both a and c are greater than 1, we use 
integration by parts at the upper vertex according to fig. ^ and obtain 

Ga,b,c = Ga-l,b+l,c — Ga,b+l,c-l- (4.3) 
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a- 1 



c- 1 



b+\ 



b+1 



(a) (b) (c) 

Figure 6: Integration by parts. 



This can be repeated until either a = or c = 

c 



c+n 



n=l 



ra=l 



a + c — n — 1 
a- 1 

a + c — n — 1 
c- 1 



G 



n,a+b+c—n,0' 



(4.4) 



The right side of this equation can be translated immediately into (" values 
by Q. For 6 = 0, we obtain the relation 



c{a,b) = {-!)' [ j;(-ir 



ra=l 



a + b — n — 1 



a - 1 



Cin)Cia + b-n) + 



n=l 



a + 6 — n — 1 
6-1 



C(n,a + 6-n) . (4.5) 



The divergent C(l) and C{l,a) appear here always in the combination 

C(l)C(a + 6-l)-C(l,a + 6-l), (4.6) 

which can be reexpressed in terms of convergent sums by the reflection 
identity (|4.1| ). In this way we arrive at 





C(a,6) = (-1)^ [ ^(-1)" 



n=2 



a + b — n — 1 



a - 1 



C{n)C{a + b-n) + 



a + 6 — n — 1 



n=2 

a + 6 - 2 
a - 1 



C{n, a + 6 — n) 



6-1 

(C(a + 6) + C(a + 6-l,l)) . (4.7) 
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Here on the left side we must assume a > 1 if regularisation is to be avoided. 
4.2 Length three 

Proceeding to sums of length three, let us again begin by exploiting the 
triviality of g^^^ . Since there are two g^^^ propagators we have now several 
possibilities. Fig. ^ shows the "standard" sea shell diagram, representing 
C(a, 6, c), as well as the diagrams related to it by a change of direction of 
one or both of the g^'^^ propagators. 



b b b b 




(a) (b) (c) (d) 



Figure 7: Diagrams related to C(a,6, c). The propagators without label are 
gio) propagators. 

Those diagrams represent the quantities 

(a) : C{a,b,c), 

(b) : C{b,a,c) + C{b,c,a) + C{b,a + c), 

(c) : C{a,c,b) + C{c,a,b) + C{a + c,b), 

(d) : C(c,6,a). 

(4.8) 

Adding those diagrams in pairs to create ^^'^-''s one obtains the following 
four identities, 

Cia,b,c)+C{b,a,c)+C{b,c,a) + C{b,a + c) + C{a + b,c)-C{a)C{b,c) = 0, 

(4.9) 
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C(a,6,c)+C(a,c,6)+C(c,a,6) + C(a + c,6) + C(a,6 + c)-C(c)C(a,6) = 0, 

(4.10) 

C{b,a,c)+C{b,c,a)+C{c,b,a) + Cib,a + c) + C{b + c,a)-Cic)C{b,a) = 0, 

(4.11) 

C{a,c,b)+C{c,a,b)+C{c,b,a) + C{a + c,b) + C{c,a + b)-C{a)C{c,b) = 0. 

(4.12) 

These identites generalize the reflection identity ( |4.5| ). Let us call them 
"permutation identities" . 

At length three we can also already make use of the three-point identities 
( p.lTp . We consider again diagram ^a, and apply the first one of the three- 
point relations to the two g^^^ propagators running into the "root" vertex. 
The result is the diagrammatic identity shown in fig. ^. 




Figure 8: Diagrammatic identity derived from the three-point identity. 



It can be translated term by term into the following identity, 

C(a,6,c) = -C{b,c,a)-Cic,a,b)+Cia + b + c) + Cia)Cib,c) 

+C(6)C(c, a) + C(c)C(a, 6) - C(a)C(&)C(c). (4.13) 
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The identities ( plO| ) through ( ^ ) can be obtained by applying the per- 
mutations 6— >a—>c—> 6, a<->c, on (|4.9| ). Taking the latter identity 
and applying all permutations, we obtain 6 identities which can be written 
as 

Mz = a, (4.14) 

where z = {({a, b, c), ({a, c,b), . . . , C(c, b, a)) (all permutations of the argu- 
ments) and a is a vector which contains only C values of length 1 and 2. 
The rank of the coefficient matrix M is 4, i.e. four zeta values in z can be 
expressed by the other two (and lower-length zeta values), e.g. by C{a,b,c) 
and C(0) c, b). Taking into account also relation ( [4. 131) the coefficient matrix 
M gets more rows but the rank does not change. Therefore this relation is 
up to lower- length identities not independent from ( [4.91) . When two of the 
arguments a, b, c coincide, say b = c, i* becomes {Cia, b, b), ({b, a, b), ({b, b, a)) 
and the rank of M reduces to 2. For all arguments coinciding the rank of 
M is 1. 

There is another possibility to derive an identity from diagram (^) which 
generalizes immediately to ( values of larger length. Every g^^^ propagator 
in (^) can be replaced by the reverted propagator according to 

g^S = -g^2i +Su-i. (4.15) 

This leads to a sum of diagrams where all ocurring g^'^^ propagators are 
directed towards the "root" vertex. The resulting identity 

C(c,6,a) = Cia,b,c)-Cia,b)C{c) + C{a,b + c) + 

+ C(a)(-C(6,c) + C(6)C(c)-C(& + c)) + 

+ C(a + b,c)- C(a + 6)C(c) + C(a + 6 + c) (4.16) 

can be obtained also by subtracting ( [4. 91 ) from (|4.12| ) and applying appro- 



priately the length-2 identity (4.1), thus it is not a new identity. 

Partial integrations yield additional identities. As for C values of length 2, 
we consider the more general diagrams fig. ^ which evaluate to the numbers 
Ga,k,b,i,c = Ga,k,b,c,i- Some of thesc numbers can be identified with C values: 

Ga,o,b,o,c = Gafl,b,cfi = C(«) b, c), Go^k,b,o,c = Go^k,b,c,o = C{k)C{b, c). (4.17) 

We start with ({a, b, c) = Gafl,b,o,c- First, we integrate by parts at the upper 
right vertex until 6 = or c = 0. The combinatorics is the same as in eq. 



(4.4). The terms with c = can be identified with C values. In the terms 



where 6 = 0, the two inner propagators can be exchanged according to the 
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Figure 9: Auxiliary diagrams for partial integrations at length 3. 




(a) (b) (c) (d) (e) 



Figure 10: Identity which proves that k and / can be exchanged when the 
upper propagator is zero. Diagram (c) is zero because no propagator goes 
into the upper right vertex. Diagram (e) is zero for the same reason at the 
right vertex. 



identity depicted in fig. 10. After this exchange we can integrate by parts 



at the upper left vertex until terms with a = or = arise. In terms of (" 
values, we obtain 



C{a,b,c) = 

= (-i)^E 



n=l 



E/a + b + c - m - n - 1 . ^. , , 

](,(m,a + + c — m — n,n) + 
V b — l,a — m,c — n 

.771=1 ^ 



b+c—n 

+ E (-1) 

m=l 



b + c — n — l\ fa + b + c — m — n — 1 



6-1 J\ a-l 

'xC,[m)C,{a + b + c — m — n,n) 



+ 



n=l 



c-1 



(4.18) 



where on the left hand side we have again to require that a > 1. The 
divergent terms for m = 1 appear always in a combination which can be 
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eliminated by identity ( 4.10| ), 

C(l)C(a, b) - C(l, a, b) = C(a, 6, 1) +C(a, 1, 6) +C(a + 1, + C(a, !)• (4.19) 

Alternatively one may after the first step, in the terms with 6 = 0, rather 
than interchanging the two inner propagators, those labelled "k" and "1" in 
fig. 9, instead interchange propagators "k" and "c" . Proceeding in the same 
way as before one obtains the following identity, 

C(a,6,c) = (-ir^('^ + ^^-^-\(a,n,6 + c-n) 

n=l ^ ^ 

n — 1 m — 1 ^ ' 



b + c — n — 1\ fa — m + n — 1 
, , n — 1 

71=1 m=l ^ ^ ^ 

x^(m, a — m + n,b — n + c) 

a — m + n— l\/6 — n + c— 1 
a- 1 /V b-1 



n=l m=l 

xC(m)C(a - m + n, 6 - n + c). (4.20) 

Here again the terms involving a ^(1, . . .) appear only in the combination 
(P^) and thus can be removed without the need for regularisation. 



4.3 Arbitrary length 

All three different procedures which we have used for constructing multiple 
C identities - partial integrations, reversion of g^^^ propagators, and the use 
of the three-point identity - can be generalized to the arbitrary length case 
without difficulty. 




Figure 11: Diagrammatic representation of permutation identities for arbi- 
trary length. The propagators without label are g^^^ propagators. 
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The generalisation of the permutation identities to an arbitrary length 
is based on finding representations of sea shell diagrams with reverted inner 
^(0) propagators, which are added in order to use the two-point relation 
( [4.15| ). Let us consider first the simplest case with one reverted propaga- 



tor (fig. O). Diagrams ll.b-d can be expressed immediately by C values. 



Diagrams |ll|.a-c represent the series 



ni>...>nm>0, 
nm+l>--->n.;>0 



(4.21) 



By decomposing the summation range into regions where the tia's are com- 
pletely ordered, we obtain a sum of C values which can be constructed as 
follows. 

Let us denote by Ma^^a2,a^2 C {{I}, {2}, {1,2}}"^+"^+"^^ the set of all 
(oi +a2 -|- ai2)-tuples consisting of ai elements {1}, 02 elements {2} and 012 
elements {1, 2}. Further, we define a map (1 = 111 + m') 

{ki , . . . , km', ki; nil 1 ■ ■ ■ : ITT'l—a) ' * 

{h+b[,...,bi.a + b'i_^), (4.22) 

where (61, ... , bi-a) results from (mi, . . . , mi-a) by replacing {2} with 0, and 
{1} and {1, 2} with ki, . . . , km (in this order); (6'^^, . . . , b[_^ results from 
(mi, . . . , rrii^a) by replacing {1} with 0, and {2} and {1, 2} with km+i, ■ ■ ■ , 
ki (in this order). For example, 

p\Jki, k2; fcg, k^; {1}, {1, 2}, {2}) = [ki, k2 + ks, k^). (4.23) 



With these definitions, the identity depicted in fig. 11 can be written as 

mm(m,i— m) 

Yl C{Pm,l-m(.h,---,km;km+l,...,kr,C)) = (4.24) 

a=0 ^&Mm,-a,l-m-a,a 

Ciki, . . . , km) C{km+1, ■ ■ ■ , k). 

The generalization to the cases with more than one reverted propagator 
should be obvious. But we note, that in these cases the C values of maxi- 
mal length appear in combinations which can be constructed also from the 
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identities ( 4.24| ). Thus we conjecture that ah identities which are based on 



the reordering of summation ranges are generated by (f4.24| ). 



Imitating the considerations on page 16, we can write all permutation 



identities ( 4.24 ) for fixed / (but varying m), where the arguments of the 



length-/ ^ values are taken from the set {ki, ... ,ki} in ail possible orderings, 



in the form (|4.14D , where now the components of z consist of all different 
values which result from C{ki, ■ ■ ■ ,ki) by permutations of the arguments, a 
contains only ( values of length less than Assuming ki, . . . , ki mutually 
different, we found that the coefficient matrix M has rank 18 for / = 4 
and rank 96 for 1 = 5. These results suggest that for length / the rank 
is /! — (/ — 1)! and that the permutation identities suffice to express the l\ 
considered C values of length / by the subset of (/ — 1)! C values where one 
of the arguments is held fixed at a certain position. 

To generalize the partial integration procedure from length three to 
length m, we can proceed in various ways. For example, we can simply 
iterate the above exchange of inner propagators. In the first step one ap- 
plies the same partial integration as in the length three case to the rightmost 
vertex of the sea shell diagram (fig. |2|). For those terms in the result where 
km-i = one makes use of this new g^^^ propagator to interchange the 
adjacent inner propagators. Then one performs partial integrations on the 
next-to-rightmost vertex until either km-2 = or the next inner propagator 
becomes the zero propagator. In the first case the procedure continues with 
another interchange of inner propagators. After maximally m — 1 such prop- 
agator interchanges, proceeding from right to left, the final step is reached, 
which is again the same as in the length three case. 

In the first step we have the same ambiguity as before. Depending on its 
resolution we arrive at a generalization of either ( 4.1S| ) or ( 4.20| ). We give 
here the formula generalizing eq.( [4.20| ). 



c{k^,...,k^) = {-If- 'y. C"''' IZ-T^ 

, i\fcm f k'm—l — Um-l -\- n^ji — \\ / krn-2 ~ "^^^,-2 ~'r Tlra-l 

nm-i=lnm-2=l 

^C(^l; • • • ) ^m— 3j T^m—2i km— 2 ~ '^m—2 "I" ^m— 1) ^m— 1 ~ T^m—l + ^m) 
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H-i)'- E E E 

"m-l=l nm-2 = l »lm-3 = l 

km-3 — nm-S + f^m~2 — 1 
nm-2 - 1 



nm~l km-3 /, I 1 \ /i- _l_ 1^ 



km—l 1 / \ km— 2 1 y 

C(A;i, . . . , km-4,'nm-3, km-3 " "-m-S + "-m-2, • • •) 



rtm 713 kl m—1 , i\ /T , IN 

+ ^ ^ E E n ' """"" ''^"""^ 

nm-i=lnm-2=l n2=l'^i=l™=2 

xC(ni, fci - ni + 722, . . . , /Cm-l - ?T-m-l + ?T-m) 
nm nm-1 n2 m—1 

+ E E •••E(-i)'^'"-"^n 

nm_i=lnm_2=l ni=l n=l 

xC{ni)C{ki - 77-1 + 722, fc2 - 722 + ns, . . . , ^m-l " nm-1 + ?^m), 



- 1 /V "2 - 1 

to) 



kn 1 



where 71^, = km- Let us also give the special case km = 1 of this formula 
which is particularly simple, 

aku...,km-i,l) = ai)Ciki,...,km-i)- (4.26) 

m—1 k^ 

~ ^ ^ ^ ^ C(^l) • • • ! k^-l, kfi -\- 1 — 72k, ^k, ^K+l, • • • , km—l)- 

K=l nre = l 

The terms involving (^(1, . . .) can be removed by means of a special case of 
eq. ( [4.24 ), namely 



C(l) Ciku • • • , km-i) - C(l, ki,..., km-i) = (4.27) 



m—1 



E [c(^i,---,^K-i, + 1, ^K+1, • • • , km—l) + 



K=l 

+ C(^l ) • • • ) ^K; 1) ^K+l , • • • , fcm— l)] • 

Up to here, the partial integrations were sequentially applied starting 
at the rightmost vertex and ending at the leftmost one. An interesting 
alternative is to do the opposite. Consider fig. ^ (page |^). As in the 
derivation of the other partial-integration identities, this diagram represents 
({ki)({k2, ■ ■ ■ , km)- On the other hand, we can use partial integrations at 
the leftmost vertex until we have only terms where either ki or /c2 became 
0. If /ci = then the resulting diagram represents a multiple C value. If 
k2 = then we exchange the adjacent inner propagators and repeat the 



(4.25) 
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whole procedure at the next-to-leftmost vertex. This continues until the 
rightmost vertex is reached, where after the partial integrations all terms 
represent multiple C values. The result is 

^2 / -\- h ft l\ 

^ k -I jC{h + ki-ni,ni,k3,...,km) + 



n\=\ 

m—1 no ni k^+i k 

X 



+ - nx-l - 1 



EEE E En 

K=2ni=ln2=l nK-i=l nft=l A=2 

Uk-i - 1 

C(^2 + - ni, ks + rii - n2, ■ ■ ■ , fc^+i + n«,_i - n^, n^,, k^,+2, ■■■,km) + 

no ni nm-2 m /, , 

ni=ln2=l nm_i=lA=2 V ^ / 

C(A;2 + no - ni,/c3 ni - ^2, . . . , nm-2 - n„_i,nm-i), (4.28) 

(no = ki). The right hand side contains no divergent terms when ki,k2 > 2. 
For /ci = 1, A;2 > 2 it becomes finite when combined with (|4.27| ). 

From our derivation clearly one would expect eqs. ( 4.25| ) and ( |4.28| ) to 



be equivalent. And indeed, it is a matter of pure combinatorics to show 
that ( |4.25| ) becomes trivially fulfilled if ( [4.28|) i s used on the right hand side. 
Nevertheless, we suspect that the form ( [4.28|) may be more useful for the 
application of these formulas to the problem of constructing a minimal basis 
of independent multiple ^ sums. 




Figure 12: Modified sea shell diagram. 
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Figure 13: "Peacock" diagram. The propagators without label are g^^^ 
propagators. 



Eq. ( ^.28| ) is a special case of a class of identities derived from (weight- 
length) shuffle algebras. In order to derive that whole class, consider Fig. 
p!3[ It evaluates to a number which we denote by 

Z{ai, . . . ,am\bi, . . . ,bn\ci, . . . ,Ck) = Z{ai, . . . , o^lci, . . . , c^l^i, . . . , 

(4.29) 

Partial integrations at the top vertex yield, similarly to ( [4.4[ ), 
Z(. . . ,am,0|6i,62, • • • |ci,C2, . . .) = 

b 



^ /&i + ci^ ^ j2'(. . . ,am,bi + ci - u\iy,b2, ■ ■ ■ |0,C2, . . .) 
+ XI ~^bi- 1 • • + ci - u\0,b2, . . . \iy,C2, . . .). 



(4.30) 



Considerations like in Fig. |10| show 

Z(ai, . . . ,Om|0, 6i, . . . |ci, . . .) = Z{ai, . . . ,am\ci, . . .\0,bi, . . .) 

= Z{ai, . . . ,ajn,0\bi, . . .\ci, . . .). 

(4.31) 
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Starting with 



C(ai, ...,am) Cih, ...,bn) = Z{0\ 



(4.32) 



and applying continuedly eqs. (|4.30D and ( 4.31| ), we end up with terms of 
the form 



Z{ai,...,am\0\bi,...,bn) 



Z{ai, . . . ,am\bi, . . . ,bn\0 
C(ai, . . . ,am,bi, . . . ,bn). 



(4.33) 



The resulting identities have the same form as the shuffle identities in the lit- 



erature, but here advantageously the binomial coefficients in ( 4.30 ) explicitly 
encode (in part) the combinatorics implicit in the shuffle algebra. 

Additional multiple C identities can be derived using the three-point 
identity ( |3.17 ). At low lengths and levels it turns out that the multiple ^ 
identities obtained in this way are not independent from the set of equations 
generated by the propagator reversions and partial integrations. Whether 
this property holds true in general we do not know. 



5 Discussion 

In the present work we have established a novel representation of multiple 
C sums in terms of Feynman diagrams in a 1 -|- dimensional quantum field 
theory. We demonstrated the usefulness of this representation for the deriva- 
tion of identities between such sums. The encoding into Feynman diagrams 
proposed here provides a very convenient book-keeping device for certain 
formal manipulations performed on such sums, as our examples should have 
amply demonstrated. 

Concerning the novelty of the identities derived here, the length- two 
identities presented in section 4.1 are, of course, well-known. At length 
three, the identities derived by the partial integration procedure, ([4. 191) , 
( ^1.201 ), are similar, and presumably equivalent, to the 'decomposition equa- 
tions' derived in by explicit series manipulations (their eq.(l)). Similarly, 
their 'permutation equations' (2) coincide with our eq. ([4.10| ). Eq.( 4.13| ) is 



contained as a special case in Theorem 2.2 of |21]. However, we have not 
been able to locate in the literature an exact equivalent of our length m 
identities (t4.25| ),( 4.28|) 0. The only identities available for arbitrary lengths 



and levels are those based on the 'shuffle algebra' |2^, ^ and its 

generalizations |23]. Of those the 'depth- length' shuffle identities (which 



3 



The special case obtained by combining eqs. (4. 26) and (4.27) is Theorem 5.1 of |2 
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are also called 'stuffle identities' or '* products' |^) are obviously related, 
and in fact equivalent to our 'permutation' identities, as we have convinced 
ourselves. Similarly the 'weight-length' shuffle identities are clearly related 
to our various 'partial integration' identities. In this case the question of 
equivalence is more difficult and will require further investigation. 

The reader will have noted that we did not make use at all of the precise 
form of the path integral action. Our considerations required the presence 
of all propagators 5^*^^ as well as of all the vertices V^''^, however they did 
not determine the statistical weights with which they should appear in the 
Feynman diagrams. Our choice of the weights for the propagators is mainly 
motivated by the fact that it leads to a suggestive form for the free path 
integral determinant. Namely, a simple application of the "Indet = tr In" 
identity shows that, formally, 

lnZ(0,A) = const. + ^A"^. (5.1) 

n=l 

Comparing this expression with the well-known formula for the logarithm 
of the r function 



r(l + x) = exp 



^ ' T7. 



n 

n=2 



(5.2) 



we see that we can identify the free partition function with the T function 
under the assumption that the ill-defined (^(1) appearing in ( |5.1[ ) is renor- 
malized to Euler's constant 7,^ 

■Z'renorm(0, A) = COnst. X T{1 - X). (5.3) 

Similarly the total propagator becomes relatively simply. Using the integral 
representation eq. ( p.llj) of the polylogarithm it is easily shown that 



Pl2 



J2 A'=(2-i)M? = gfS + ^ 2Fi(l, 1 - A; 2 - A; z,,). 
k=0 



(5.4) 



^This calculation may be seen as a "chiral" generalization of the calculation of the 
Scalar QED Euler-Heisenberg Lagrangian performed in |3q ]. 

^Considering the identities 7 = lim„—,cxj (X]fc=i i ~ Inri) and X^fe^i p" ~ + 7 + 0{e) 
this assumption seems quite natural. 
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For the interaction term there seems to be no such preferred choice. A 
question of obvious interest (but equahy obvious difficulty) is whether non- 
trivial interaction potentials V{g, g) exist such that the C-model would be 
exactly solvable. 

Recently the following generalization of Euler-Zagier sums (^]^) has at- 
tracted some attention, 

ni>n2> - >nm>0 "l ""^ 

where aj = ±1. Those "alternating" Euler-Zagier sums arise naturally in 
the calculation of ultraviolet divergences in renormalizable quantum field 
theories, and in the application of knot theory to the classification of those 



divergences [13]. At the same time the inclusion of alternating Euler-Zagier 
sums seems to simplify the problem of reducing the set of all such sums to 
a basic set via multiple C identities [10|. (For a tabulation of alternating 



series see [gj.j 

More generally, arbitrary A^-th roots of unity in place of the aj have been 
considered in connection with the study of mixed Tate motives over Spec 



Z 1 18, 33|. Those phase factors can be easily accomodated in the C-model. 
To generate a phase factor a = e^'^'* the propagator (|2.1[) has to be simply 
replaced by 

°° p27ri(n+s)ni2 
n=l (2^1^) 

On the path integral level this can be achieved by changing from periodic 
to twisted boundary conditions, x{l) = ax{Q), and replacing d hy d — 27ris. 
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